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Background: Boundary Layer Ingestion 


• Smith, L. H., “Wake Ingestion Propulsion Benefits,” AIAA Journal 
of Propulsion and Power, Vol. 9, No. 1, Jan-Feb 1993, pp. 74-82. 



® Boundary Layer Ingestion 
(BLI) Propulsion has potential 
for significant reduction 
(5%-1 0%) in aircraft fuel burn. 


• BLI may present significant flow distortion to fan. Aeromechanical 
response of fan in presence of distorted flow must be understood, 
and applied design must be aeromechanically robust. 
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CFD & Aeromechanics Code: TURBO-AE 


• Unsteady Reynolds-Averaged Navier Stokes equations 
® Characteristics-based finite-volume solver 
® Newton-iterative, implicit, time-accurate 
a Structured multi-block code 
a Decoupled k - e turbulence model 
a Sliding interfaces where applicable 
a Inlet distortion boundary condition 
a Throttle exit boundary condition 
a Dynamic grid deformation for blade vibration 

• One-way prescribed harmonic biade vibration (previously validated 
in aeromechanics analyses in clean flow and applied herein) 

• Two-way coupled blade motion (implemented and applied herein) 
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Computational Research Fan 





Passage Mesh 214/76/58 

Blade Mesh 81/70/58 

Inlet Duct Length^ 1 .387 

Exit Duct Length^ 2.365 

^ axial chords at tip radius 
^Part-speed condition studied here 


Tip Diameter 

22 in. 

Blade Count 

22 

Rotor Speed^ 

10800 rpm 

Blade Natural Frequency 

220 Hz 

Rotational Frequency 

180 Hz 

Time Steps/Passage 

15 

Time Steps/Revolution 

330 

Time Steps/Oscillation 

270 

Subiterations, Clean 

6 

Subiterations, Distorted 

9 

Clock per step. Clean 

1m:16s 

Clock per step. Distorted 

1m :53s 

Clock per rev, Clean 

6h:58m 

Clock per rev, Distorted 

10h:21.5m 
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Inlet Total Pressure 





• Clean Inlet ® Sinusoidal Inlet Distortion 

(4% mean-to-peak) 

Azimuthal orientation: 0° at 12 o’clock, increasing clockwise. Fan 
rotates counterclockwise. 
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Rotor Part-Speed Performance Characteristics 
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® Speedlines for clean and distorted inlets. 

» Averaged over one revolution. 

® Distorted flow yields slight degradation in mass flow and pressure 
recovery. 
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Blade vibrations are prescribed in a selected mode frequency {f), 
and nodal diameter pattern (ND) or inter-blade phase angle (a). The 
work done on the vibrating blade is calculated as follows: 


where p is the aerodynamic pressure, A is the blade surface area 
vector, X is the displacement vector on the blade surface, and t 
denotes time. The aerodynamic damping ratio (C) can be 
approximately related to the work-per-cycle ( W) and the average 
kinetic energy {Ke) of the blade over one cycle of vibration through the 
following expression: 



pdA- — dt 
ot 


( 1 ) 


surface 



( 2 ) 
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One-way Coupling: Time Marching Scheme 



Aerodynamic damping is calculated from the aerodynamic work 
quantity. Positive aerodynamic work signifies negative aerodynamic 
damping and indicates the possibility of flutter. 
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Normalized Aerod>matnic Work [-] 
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One-way Coupling: Clean Inlet, Near Op-Line 


Convergence of Aerodynamic Work for All 22 Blades 
(T = 0° for 0 ND, and a = 180° for 1 1 ND 



Oscillal.ifm Cycle Cotml Azimulhal Position of Root L^E at end of oscillation cycle [cleg] 

® May require a few cycles to converge (1 1 ND here) 

® Consistent about annulus (clean inlet) 
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One-way Coupling: Distorted Inlet, Near Op-Line 


Convergence of Aerodynamic Work for All 22 Blades 
(7 = 0° for 0 ND, and a — 180° for 1 1 ND 




® May require a few cycles to converge (1 1 ND here) 

9 Work is a function of azimuthal position: 

• Unique for each subset of blades within each ND pattern 

• May exhibit mixed stability/instability (0 ND here) 
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One-way Coupling: Aerodynamic Work Versus 
Azimuthal Position 



Azimuthal Position of Root LB at end of oscillation cycle |dcg| Azimuthal Position of Root LB at end of oscillation cj'cle |dcgj 


cT = 16.36° 

1 ND 

avg blade, clean 

-0.068% 

all blades, clean 

mixed 

avg blade, distorted 

-0.072% 

all blades, distorted 

mixed 


(T = 32.73° 

2 ND 

avg blade, clean 

+0.49% 

all blades, clean 

mixed 

avg blade, distorted 

+0.09% 

all blades, distorted 

mixed 
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One-way Coupling: Nodal Diameter Sweep 



• Aerodynamic damping versus nodal diameter using average-blade 
aerodynamic work from one-way coupled method. 
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Two-way Coupling: Mathematics 


From the original equation of motion, 

mm + ic\m + [K]{x} = {F} - {/,} (3) 

where 

{F} = - j pdA (4) 

surface 

Pre-multiply by the mode shape [O]^, establish the generalized 
displacement coordinate rj, and normalize to unit modal mass: 

[0]’'[M][<l>] = 1 (5) 

{X} = |0]{,} (6) 

Applying the conservative simplification of neglecting material and 
structural damping (i.e., assuming C = 0) yields independent second 
order ordinary differential equations for each mode /: 

Vi + ^fvi = [^/] ^{Fi- fs,i} (7) 
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Two-way Coupling: “Static Modal Force” 


• Static modal force: Modal force on rigid blades due to unsteady 
pressure loading of clean and distorted flows. One blade shown. 

= - / •^J pdA ( 8 ) 

surface 
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Two-way Coupling: Time Marching Scheme 



Aerodynamic damping is quantified by logarithmic-decrement analysis 
of blade modal displacements. Increasing magnitudes of 
displacements signify negative aerodynamic damping and indicate 
possibility of flutter. 
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Two-way Coupling: Near op-line, 0 Nodal Diameter 


All 22 Blades, a = 0° 


• Clean inlet 


• Distorted inlet 




® Decreasing magnitudes of oscillation indicate flutter stability. 

® Damping: One-way = +1 .64%; Two-way = +1 .60% to +1 .80%. 
® Clean and distorted inlets nearly-identically stable. 
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Two-way Coupling: Near stall, 1 ND FTW 


• Clean Inlet 


All 22 Blades, a = 16.36° 

® Distorted Inlet 



O.OOC8 


0.0(K16 


S 0.tX3(M 

g 

L> 

OiX ) C 2 

h 

"O 

tu 


O - 0.0002 


- 0,0004 


- 0.0006 


- 0.0008 



Time [s\ 


0.0008 


omm 


0.0004 


{) um2 
a 


o -0.0002 


- 0.0004 


- 0.0006 

- 0.0008 



Time [sj 


One-way avg blade 

-0.068% 

One-way all blades 

mixed 

Two-way coupled 

Si -0.038% 


One-way avg blade 

-0.070% 

One-way all blades 

mixed 

Two-way coupled 

Si -0.20% 


National Aeronautics and Space Administration — www.nasa.gov Fan Flutter / Inlet Distortion / Two Coupled Methods 17/23 




50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference 


Cleveland, Ohio 


Two-way Coupling: Near stall, 2ND FTW 


All 22 Blades, a = 32.73° 


® Clean Inlet 


® Distorted Inlet 
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Summary 




• A two-way coupled flutter analysis method has been implemented 
in TURBO-AE. 

® One-way and two-way coupled flutter analysis methods have been 
applied to study fan flutter in presence of inlet distortion (total 
pressure). 

• Distortion defined as single, once-per-revolution distribution about 
annulus with a 4% (mean-to-peak) variation about clean inlet’s 
radially-varying total pressure profile. 

« Fan run at part-speed with artificial mode shape and natural 
frequency. 
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Conclusions 


® Distorted flow conditions consistently yield slight degradations in 
mass flow, pressure recovery, and flutter stability versus the clean 
flow conditions of identical circumferential-mean total pressure 
distribution. 

• Flow nonuniformity and nodal diameter pattern determine each 
blade’s aerodynamic work response. Individual blade 
aerodynamic work response will vary consistently about the 
annulus from cycle to cycle due to interplay of blade natural 
frequency and rotor rotational frequency. 
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Conclusions - 

■ Continued 



• Two-way coupled method corroborated one-way coupled method’s 
average-blade flutter assessment when: 

• a// blades showed stability about the whole annulus using one-way 
coupled method; 

OR 

• “average” blade showed /nstability. 

• For 2 ND FTW near stall: 

• One-way coupled method’s average-blade showed stability in both 
clean and distorted flows. 

• Two-way couple method showed stability in clean flow but /nstability 
in distorted flow. 
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Future Work 




® Refined time-stepping, longer simulation times, refined gridding: 

• More clearly establish existence and magnitude of instabilities; 

• Interrogate higher-fidelity results to better understand fundamental 
aeromechanics, (incidence angle, shock position, e.g.) 

9 Improved fluid-structure time-coupling 

• More realistic, more complete inlet distortion definition (swirl, e.g.) 

• Physical experiment 


National Aeronautics and Space Administration — www.nasa.gov Fan Flutter / Inlet Distortion / Two Coupled Methods 22 / 23 




50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference 


Cleveland, Ohio 




niet Distortion / Two Coup 




